On the equation of state of a dense columnar liquid crystal 
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An accurate description of a columnar liquid crystal of hard disks at high packing fractions is 
presented using an improved free-volume theory. It is shown that the orientational entropy of 
the disks in the one-dimensional fluid direction leads to a different high-density scaling pressure 
compared to the prediction from traditional cell theory. Excellent quantitative agreement is found 
with recent Monte-Carlo simulation results for various thermodynamic and structural properties of 
the columnar state. 
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The lyotropic columnar liquid crystal state, character- 
ized by a two-dimensional hexagonal stacking of columns 
each with a liquid-like internal structure, has received 
considerable attention in recent years both in experimen- 
tal colloid science 0, IjJ and computer simulations 0, Q . 
The recently developed colloidal systems of polymer- 
grafted gibbsite platelets are known to show a first order 
phase transition from a spatially homogeneous nematic 
to an inhomogeneous columnar phase upon densification 
The columnar signature of the latter is evidenced 
by its bright Bragg-reflections for visible light and it has 
been confirmed on a more rigorous basis using Small- 
Angle X-ray Scattering (SAXS) measurements Col- 
loidal platelets constitute an emerging field of interest 
and future efforts can be put into reducing their polydis- 
persity and manipulating columnar textures by means 
of a magnetic field, both aimed at making high-quality 
single-domain columnar structures. These may be candi- 
dates for the production of e.g. colloidal photonic crystals 

i 

In this Letter we consider a simple but accurate the- 
oretical description inspired by cell theory, which was 
first applied to spatially ordered liquid crystals by Tay- 
lor, Hentschke and Herzfeld ,6„i7,]. By incorporating the 
rotational freedom of the disks into the framework we are 
able to quantitatively account for both thermodynamic 
and structural properties of the phase along its entire sta- 
bility range. The theory is shown to be universal, i.e. in- 
dependent of the generic shape or the length-to-diameter 
ratio of the disks. For this reason, the results are ex- 
pected to be significant for ongoing experimental studies 
on colloidal platelet systems. In particular, knowledge 
of the osmotic pressure may be useful for analyzing sedi- 
mentation experiments 's'l whereas structural predictions 
may prove helpful for the interpretation of SAXS results 
on dense systems of platelets 9] . 

To describe the properties in the two positionally or- 
dered dimensions of the columnar liquid crystal we as- 
sume the particles to be confined within discrete com- 



partments represented by hexagonal tubular cells which 
form a close-packed structure according to the classical 
free- volume ("cell") theory 13, The disks are al- 

lowed to take any position within the cell but they may 
not share the cell with another particle or penetrate an 
adjacent one. The cell approach will be combined with an 
appropriate description of the one-dimensional fluid be- 
haviour along the columnar direction. We start from the 
traditional Tonks model applied to a linear fluid of N 
hard disks with diameter D and thickness L whose cen- 
tres of mass can move freely on a line with length Since 
the disks are allowed to rotate freely around their centres 
of mass, the effective excluded thickness Li j between two 
adjacent disks i and j is an orientation-dependent quan- 
tity, i.e. Li.j(r2i, 51j) > L, in terms of the solid angle 
We assign to the position of particle k on the line and 
fix the first and last particle at xi = and xjq = £, re- 
spectively. The configurational integral for this system in 
the macroscopic limit L/£ is then formally written 
as 
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with V the thermal volume pertaining to the translational 
and orientational kinetic degrees of freedom. The brack- 
ets denote an orientational average according to some 
unknown orientation distribution function (ODF) /(fi) 
which is normalized according to / f{n)dil = 1. Note 
that Qn is proportional to an A'^-dimensional free vol- 
ume with Ltot being the total occupied length for a given 
orientational configuration, expressed in terms of the fol- 
lowing sum 

N 

Ltot(f^l, ■ • ■ , ^n) = ^ Lk,k+l{^k, ^k+l)- (2) 
k=l 

Eq. ^ is difficult to analyze rigorously so we shall ap- 
proximate it as follows 

i'^ ( \ l~ \ ^ ^ 

which is assumed to be justified for the strongly aligned 
orientational configurations we expect in a dense colum- 
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FIG. 1: (a) Equations of state for the columnar phase. The sohd curve is the prediction from the present cell model, the 
dashed one is the traditional cell pressure, given by Eq. (1161 1. The symbols correspond to the simulation data from Zhang et 
al. 3 for D/L = 10 (crosses) and D/L = 20 (triangles). The inset shows the pressure in the dilute regime near the columnar- 
nematic ("melting") transition at 0* « 0.48. (b) Dimensionless chemical potential jSjl = 13^ — \n[vo/V(j)c-p\ from simulations 
(dotted curves) and cell theory (solid curves). The prediction from the traditional cell model, indicated by the dashed curve, 
is independent of D/L. 



nar state. Further simplification in this respect can be 
achieved by neglecting the dependency of the excluded 
length on the azimuthal angle of the particles. To this 
end we shall consider an effective disk thickness, which is 
determined solely by the polar deflection angle 9 between 
the particles' symmetry axis and the line unit vector. For 
small angles this quantity is given by 



LcS — L 
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up to leading order in 6. The factor '1/2' in Eq. is in- 
cluded explicitly to approximately restore the azimuthal 
dependency of L^ff. Setting the factor equal to unity 
yields the maximum L^ff at fixed polar angles, which is 
only realized if Aip for two adjacent disks is exactly tt. 
In reality, Aip is randomly distributed between and 27r 
so that the total excluded length is usually way below 
its maximum value. The orientationally averaged total 
occupied length is approximated by the following mean- 
field expression: 
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where the ODF obeys common uniaxial symmetry and 
depends only on the polar angle 9. The configurational 
integral then simply becomes 
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in terms of the linear density p = NL/£. Following On- 
sager ^13] the orientational configurational integral 
reads 
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and the total Helmholtz free energy /3-Ffluid = — In Qn 
of the modified Tonks fiuid is given by a superposition 
of the ideal, orientational and configurational entropic 
contributions: 
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with P — l/fc^T- The thermodynamic equilibrium ODF 
can be obtained by a functional minimization of the free 
energy with respect to the orientational distribution un- 
der the normalization restriction. After some algebra we 
arrive at the following closed expression for the normal- 
ized ODF: 



f(0)^—exp[-a\9\l 



(8) 



where a depends on the density and the disk aspect ratio 
D/L via 
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Since a 3> 1 for sufficiently anisometric disks {D/L ^ 1) 
the ODF is sharply peaked around 6* = 0, as we expect. 
The divergence of a at close packing (p = 1) indicates 
that the disks are forced to orient parallel to the direc- 
tor in this limit. Inserting Eq. ^ and straightforward 
integration yields an explicit free energy in terms of the 
density p. Taking the standard derivative with respect to 
the density then gives the (dimensionless) pressure pPL 
of the one-dimensional fluid: 



pPL = 3p/(l - p), 



(10) 
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which is exactly three times the original Tonks pressure 
[T3 |. This result implies that the orientational confine- 
ment of the disks inside the columns gives rise to an addi- 
tional entropic contribution Ipj (1 — p) to the total pres- 
sure. Note that the derivation of Eq. (|10l) closely resem- 
bles the origin of the quasi-ideal pressure of the nematic 
state /3P ~ SN/V within the Gaussian scaling Ansatz for 
the ODF I14lll5j . although both have a different physical 
basis. 

The description of the columnar phase in the two posi- 
tionally ordered dimensions is analogous to that of a 2-D 
(hexagonally) ordered configuration of N disks. Accord- 
ing to the non-correlated version of the cell model the 
configurational in tegr al of the A^-particle system can be 
approximated by 



Q5^" = J dr^exp[/3C/(r^)] ~ (^J dvexp 
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where it""jj(r) is the potential energy between the parti- 
cle and its nearest neighbours. For hard interactions the 
second phase space integral is simply the (2-D) free vol- 
ume of the particle in the cell. Assuming that the nearest 
neighbours constitute a perfect hexagonal cage, this free 
volume is given by Vfrce = y3(Ac - Df/2 with Ac the 
nearest neighbour distance. The configurational integral 
then becomes (ignoring all irrelevant contributions) 
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in terms of the dimcnsionlcss spacing Ac = Ac/D. Ap- 
plying the condition of single-occupancy (i.e. one array 
of disks per column) we can use Ac to relate the linear 
density p to the three-dimensional volume fraction via 



(12) 



with (h* = 



cp 



the volume fraction normalized to its 



close-packing value 0cp — tt/2^/3. 

The total excess Helmholtz free energy of the columnar 
state is obtained by adding the cell contribution to the 
Tonks excess free energy. Omitting all constant terms we 
arrive at 
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(13) 

representing the orientational, 'Tonks' and cell contri- 
butions, respectively. Inserting Eq. (|12|l and minimiz- 
ing the free energy l20| with respect to the cell spacing 
Ac yields a third-order consistency equation with com- 
plicated solutions. Expanding the physical solution near 
close-packing in terms of 1 — 0* up to leading order yields 

Ac = l + i(l-<^*) + O[(l-0*)']. (14) 
5 

Substituting its closed form into the free energy Eq. 
(|13|l and taking the appropriate derivative gives a sim- 
ilar expansion for the total dimensionless pressure P — 
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FIG. 2: Normalized inter- and intracolumnar spacings, Ac/D 
and A„/L respectively, as a function of (p* . Solid lines are 
theoretical predictions, the symbols follow from simulations 
for D/L = 10 (crosses) and D/L = 20 (triangles). The dotted 
curve follows from the traditional cell model and denotes both 
spacings. Inset: Average polar angle relative to the internal 
angle (Sint ~ L/D) plotted versus 4>* . 



(3Pvo/(j)cp (with vq the disk volume), i.e. 
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+ 6.4 + 1.128(1 - (j)*) + 0[(1 - (j)*yi (15) 



indicating that the high-density scaling pressure, i.e. the 
leading order contribution, is essentially different from 
the classical cell prediction 3/(1 — </>*) for hard spheres 
[TtL fl^ . The latter result is completely analogous to our 
result for perfectly aligned disks and can be reproduced 
directly from Eq. (|13() by omitting the orientational con- 
tribution given by the first term. The "traditional" cell 
pressure is given by 
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In Fig. 1(a) we have plotted the abovementioned cell 
equations of state, i.e. Eq. H16|l and the closed-form ana- 
log of Eq. (|15|) , along with the ones obtained from Monte- 
Carlo simulations on cut-spheres by Zhang et al. 
The quantitative agreement between the present cell de- 
scription and the simulations significantly improves upon 
densification. Above (ff « 0.8 our prediction is found to 
match the simulations within 1%. Obviously, the dis- 
crepancy is much larger in the dilute regime {4>* < 0.6) 
where the cell model, at least the simplest version con- 
sidered here, is no longer quantitatively reliable. 

The chemical potentials are depicted in Fig. 1(b). 
These are readily obtained from the absolute free energy 
and pressure of the cell model via the Legendre trans- 
form p — {F + PV) /N and by means of a Gibbs-Duhcm 
integration of the pressure fits obtained from the sim- 
ulation data Q. Again, close agreement is found be- 
tween both. The values at melting are found to differ 
only by Ap w 0.3fcBT and OMbT for D/L = 20 and 
D/L = 10 respectively, indicating that the cell predic- 
tion is surprisingly efficacious throughout the 0*-range. 
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The "splitting" of the curves for both aspect ratios is 
simply due to the orientational entropy [second term in 
Eq. (O], PFor/N = 2 In a - 2, which depends expHcitly 
on D /L via Eq. © and therefore gives rise to a different 
intercept at the melting volume fraction. 

Let us now focus on the spacings between the columns 
Ac and the average disk spacing inside the columns A„ = 
A„/L = p~^, which follows directly from Eqs. H12|) and 
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Comparing with Eq. (|14() we see that the intracolum- 
nar distance A„ between the disks grows faster than 
the intercolumnar one Ac. The expansion of the colum- 
nar structure is therefore anisometric. This behaviour 
is quite different from the classical model for which 
Ac = A„ = (f)*~^/'^ indicating an isometric expansion 
upon lowering <j)* . In Fig. 2 the predicted spacings, given 
by the closed- form analogs of Eqs. (|14ll and (|17ll . are 
compared with simulation results. Unlike the pressure, 
the prediction for the columnar spacings remains surpris- 
ingly accurate even in the regime close to the columnar- 
nematic transition. Quantitatively, theory and simula- 
tion are found to agree within 1% over the entire colum- 
nar stability range. 

To assess the degree of orientational order in the dense 
columnar state we consider the ratio of the average po- 
lar angle {\9\) to the internal angle O-mt ~ L/D, i.e. 
(|6»|) /^int = 2D/La. Substituting Eqs. © and ^ 
reveals a simple proportionality with respect to the in- 
tracolumnar spacing, i.e. {\9\) /6'int — (4/3)(A„ — 1). 
This tells us that the disks are only marginally perturbed 
away from their parallel orientations since the average 
"off-parallel" deflection angle does not exceed the inter- 
nal angle of the disk. The full result, included in Fig. 
2, explicitly shows that this situation remains up to the 
columnar-nematic transition, located around 6* ~ 0.48. 



From this we conclude that the orientational freedom of 
the disks is extremely small throughout the entire colum- 
nar stability regime. Moreover, the dominance of near- 
parallel configurations gives an a posteriori justification 
for the mean-field asymptotic analysis presented here. 

In conclusion, we have constructed a modified cell the- 
ory for the columnar state by explicitly accounting for 
the rotational freedom of the disks. Our approach con- 
stitutes a significant quantitative improvement over the 
traditional one which, contrary to the case of a hard- 
sphere FCC-crystal 18j, appears to be inappropriate for 
a columnar liquid crystal even in the regime near close- 
packing. The theory does not only quantitatively pre- 
dict the thermodynamics of the columnar state, as evi- 
denced by the pressure and chemical potential, but also 
structural features in terms of the characteristic colum- 
nar spacings. 

Future work could be aimed at applying the present ap- 
proach to a dense smectic-A phase of hard spherocylin- 
ders by considering a two-dimensional fluid of rotating 
rods . Although the theory is expected to be less suc- 
cessful here, it would be intriguing to verify its quantita- 
tive merits for a lyotropic smectic phase. Furthermore, 
the cell description could be refined in the regime close 
to the melting transition by introducing more advanced 
cell theories which include e.g. cooperative motion of the 
columns [T9l |. The fact that the simulation pressures in 
Fig. 1(a) remain insensitive to the aspect ratio through- 
out the entire columnar stability range is surprising and 
supports the notion that the columnar phase is domi- 
nated by cell-behaviour up to the melting transition, al- 
beit in a more sophisticated fashion than we described 
here. 
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